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Abstract. Exploiting recent regularity estimates for the Monge- Ampere equation, under some suit- 
able assumptions on the initial data we prove global-in-time existence of Eulerian distributional solu- 
tions to the semigeostrophic equations in 3-dimensional convex domains . 



1. Introduction 

A simplified model for the motion of large scale atmospheric /oceanic flows inside a domain ft C M 3 
is given by the semigeostrophic equations. 

Let ft C M 3 be bounded open set with Lipschitz boundary. Then the semigeostrophic equations 
inside ft are: 



(1.1) 



' d t u g t (x) + (u t (x) • V)uf(x) + Vpt(cc) 

d t m t (x) + (u t (x) ■ V)m t (x) = 

u\ (x) = JVpt(x) 

V • u t (x) = 

u t (x) ■ fu(x) = 
^p (x) =p°(x) 



-Ju t (x) + m t (x)e3 



(x,t) £!lx (0,oo) 
(x,t) £ ft x (0,oo) 
(x,t) e ft x [0,oo) 
(x,t) e ft x [0,oo) 
(x,t) G <9ft x [0, oo) 

x e ft. 



Here p° is the initial condition for vq is the unit outward normal to 5ft, = (0, 0, 1) T is the third 
vector of the canonical basis in R 3 , J is the matrix given by 



J := 



and the functions ut, pt, and mt represent respectively the velocity, the pressure and the density of 
the atmosphere, while nf is the so-called semi-geostrophic wind0 Clearly the pressure is defined up 
to a (time-dependent) additive constant. 

JVpt an d introducing the function 




Substituting the relation itf 



(1.2) 



1 



P t {x) :=Pt(x) + -(z? + xl), 



As it will be clear from the discussion later, we do not need to specify any initial condition for u and m. 
! We are using the notation f t to denote the function fit,-). 
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the system (jl.ip can be rewritten in Q x [0, oo) as 

'd t VP t (x) + V 2 P t (x)u t (x) = J(VP t (x) - x) 
, \ V • u t (x) = 

u t [x) ■ vn(x) = 

P Q {x) = p°{x) + \(x 2 +x 2 ). 

Notice that, given a solution (P,u) of ()1.3|) . one easily recovers a solution of (jl.ip : indeed p< can 
be obtained from Pt through (jl.2p and the density mt is given by nit = d%Pt (in particular, the third 
component of the first equation in (jl.3p tells us that dtmt + (ut • V)m( = is satisfied). 

Energetic considerations (see |11[ Section 3.2]) show that it is natural to assume that the function 
Pt is convex on f2. This condition, first introduced by Cullen and Purser, is related in [144 121] to a 
physical stability required for the semigeostrophic approximation to be appropriate. If we denote with 
J2fn the (normalized) Lebesgue measure on 17, then formally pt := (S7Pt)^n (see, for example, pjj 
Appendix A]) satisfies the following dual problem 



(1.4) 



f d tPt + v • (u tPt ) = o 

U t {x) = J(x - VP t *(x)) 

Pt = (VP t ) tt JSfn 
P Q {x)=p°{x) + \{xl + x 2 



2> 



Here P* is the convex conjugate of Pt, namely 



P t *(y):=sup{yx-Pt{x)) Vy E M 

and (VP()jj«Sffj is the push-forward of the measure ££$i through the map VPt : O — >• M 3 defined as 
[(VP t ) fl J2h](A) = j2h((VP t )-V)) ^ all id 3 Borel. 

The dual problem is pretty well understood, and admits a solution obtained via time discretization 
(see [5j[l3]). Moreover, at least formally, given a solution Pt of the dual problem (jl.4p and setting 

(1.5) u t {x) := [d t VP*}{VP t {x)) + [V 2 P t *](VPt(x))J(VP t (x) - x), 

the couple (Pt,ut) solves the semi-geostrophic problem (jl.ip . However, because of the low regularity 
of the function Pt the previous velocity field may a priori not be well defined, and this creates serious 
difficulties for recovering a "real solution" from a "dual solution". 

Still, a recent regularity result |15j can be applied to show that the map Pt is W 2 ' 1 in space, so that 
we can give a meaning to the second term in the definition of ut- More precisely, in [15] it is shown 
that \D 2 u\ log^_ \D 2 u\ G L\ oc for any k, and following ideas developed in [HEU], we will be able to show 
that the function Pt is regular enough also in time, so that the couple (Pt,ut) is a true distributional 
solution of (jl.lj) . 

Let us point out that the regularity result in [15] has been recently extended, independently in [17] 
and [23], to \D 2 u\ £ L 7 , where 7 > 1 depends on the local L°° norm of logp^. However, as we will 
better explain in Remark 13.61 in our situation there is no advantage in using this improvement, since 
the fact that 7 depends on || log pt Hoopoe makes the estimates less readable. For this reason we will 
rely only on the L log L integrability given by [15] , as we previously did in pQ . 

The first existence result about distributional solutions to the semigeostrophic equation is presented 
in p], where the analysis is carried out on the 2-dimensional torus (see also [21] where a short time 



GLOBAL EXISTENCE FOR THE SEMIGEOSTROPHIC EQUATIONS 



3 



existence result of smooth solutions is proved in dual variables, and because of smoothness the existence 
can be easily transferred to the initial variables). 

The 3-dimensional case on the whole space R 3 , which is more physically relevant, presents addi- 
tional difficulties. First, the equation (11. lj) is much less symmetric compared to its 2-dimensional 
counterpart, because the action of Coriolis force Jut regards only the first and the second space com- 
ponents. Moreover, even considering regular initial data and velocities, regularity results require a 
finer regularization scheme, due to the non-compactness of the ambient space. 

Our proofs are also based on some additional hypotheses on the decay of the probability measure 
Po = (VPo)tl-^h- This decay condition happens to be stable in time on solutions of the dual equation 
(|1.4p . and allows us to perform a regularization scheme. 

It would be extremely interesting to consider compactly supported initial data po = (VPo)tt=%- 
However the nontrivial evolution of the support of the solution pt under (jl.4p prevents us to apply the 
results in [15] (which actually would be false in this situation), so at the moment this case seems to 
require completely new ideas and ingredients. 

Definition 1.1. Let P : x [0, oo) — > M. and u : Q x [0,oo) — > M. 3 . We say that (P, u) is a weak 
Eulerian solution of (jl.3p if: 

- |u| G Lj£.((0,oo),L£ }C (n)), P G L£ c ((0, oo), W^°{Q)), and P t (x) is convex for any t > 0; 

- For every <fi G C™(Q x [0, oo)), it holds 

(L6) oo 

J J VP t {x)\d t (t)t{x)+u t {x)-V(t)t{x)^ + j\vPt{x)-x^(t)t{x)dxdt + J VP (x)cP (x)dx = 0; 

- For a.e. t G (0, oo) it holds 

(1.7) [ Vip(x) ■ u t (x) dx = for all ip G c7 c °°(0). 

Jn 

Remark 1.2. This definition is the classical notion of distributional solution for (jl.3p except for the 
fact that the boundary condition UfVn = is not taken into account. In this sense it may look natural 
to consider ip G C°°(f2) in (jl.7p . but since we are only able to prove that the velocity ut is locally in 
L , Equation (|1.7p makes sense only with compactly supported ip. On the other hand, as we shall 
explain in Remark 1 1A\ we will be able to prove that there exists a measure preserving Lagrangian flow 
Ft : £1 — > £1 associated to Ut, and such existence result can be interpreted as a very weak formulation 
of the constraint UfVQ = 0. 

As pointed out to us by Cullen, this weak boundary condition is actually very natural: indeed, 
the classical boundary condition would prevent the formation of "frontal singularities" (which are 
physically expected to occur), i.e. the fluid initially at the boundary would not be able to move into 
the interior of the fluid, while this is allowed by our weak version of the boundary condition. 

We can now state our main result. 

Theorem 1.3. Let ft C BL 3 be a convex bounded open set, and let J£q be the normalized Lebesgue 
measure restricted to O, that is JHn(Q) = 1. Let po be a probability density on M 3 such that po G 
L°°(M 3 ), 1/po G P^ C (M 3 ) and 

limsup (po(^)l^l^) < co 
\x\— >oo 
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for some K > 4. Let p t be a solution of f 1 1 . 4 [> given by Theorem \3.1\ P t * : ]R 3 — > R the unique convex 
function such that 

P t *(0) = and (VP;) fl ( A Jz? 3 ) = if n , 

and lei Pj : M 3 — >• M be its convex conjugate. 

Then the vector field ut in (jl.5p is well defined, and the couple (Pt,ut) is a weak Eulerian solution 
of (jl.3p in the sense of Definition \l.l\ 



Remark 1.4. Following Cullen and Feldman one can give also a notion of Lagrangian solution of 
the semigeostrophic equation. More precisely they show the existence of a measure preserving flow 
Ft : ft — > ft which solves a sort of Lagrangian version of (jl.ip (see [12] and [H Section 5] for a more 
precise discussion). Actually the flow they constructed has the explicit expression F t = VP^oGjoVPo, 
where Gt is the regular Lagrangian flow associated to the BV vector field Ut = J(x — VP/), in the sense 
of Ambrosio, Di Perna and Lions (see [3J SI [18] ) . In [TJ Section 5] we showed, in the two dimensional 
periodic setting, that for almost every x the map t \— > Ft(x) is absolutely continuous with derivative 
given by ut(Ft(x)). The proof of this fact can be almost verbatim extended to our contest, showing 
that, for almost every x £ ft, 1 1— > Ft{x) is locally absolutely continuous in [0, oo) with derivative given 
by ut(Ft(x)). We leave the proof of this fact to the interested reader. Finally we remark that the 
uniqueness of such a flow (both according to the definition given in [12] or in [1]) is unknown. 

Acknowledgement. L.A., G.D.P., and A.F. acknowledge the support of the ERC ADG GeMeThNES. 
A.F. was also supported by the NSF Grant DMS-0969962. M.C. and G.D.P. also want to acknowledge 
the hospitality of the University of Texas at Austin, where part of this work has been done. 

2. Regularity of optimal transport maps between convex sets of M 3 

Throughout this paper, S] C R 3 is a bounded convex open set, dn > is fixed in such a way that 
ft C P(0, dn), and Jz?q denotes the normalized Lebesgue measure restricted to ft. 

In this section we recall some regularity results for optimal transport maps in ]R 3 needed in the 
paper. 

Theorem 2.1 (Space regularity of optimal maps between convex sets). Let fto, ^i be open sets o/M 3 , 
with fti bounded and convex. Let p = pJz? 3 and v = cr=Sf 3 be probability densities such that //(f2o) = 1> 
z/(S7i) = 1. Assume that the density p is locally bounded both from above and from below in namely 
that for every compact set K C ft there exist Xq = Xo(K) and Ao = Aq(K) satisfying 

< A < p(x) < A V x e K. 

Futhermore, suppose that Ai < a(x) < Ai in Then the following properties hold true. 

(i) There exists a unique optimal transport map between p and v, namely a unique (up to an 
additive constant) convex function P* : f?o — >• M such that (VP*)j// = v. Moreover P* is a 
strictly convex Alexandrov solution of 

p(x) 



det V z P*{x) = f(x), with f{x) 



a(VP*(x)) 



P* G Wf£(n ) n (7j,f (fi ). More precisely, if ft <e ft 

is an open set and < A < p(x) < 
A < oo in ft, then for any k € N there exist constants C\ = C\{k, ft, fti, A, A, Ai, Ai), f3 = 
/3(A, A, Ai, Ai), and C2 = C2(ft, fti, A, A, Ai, Ai) such that 



[ |V 2 P*|log^|V 2 P*|dx < d, 
Jn 
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and 

\\ p *\\c l >f>(a) ^ ^2- 

(iii) Let us also assume that Qq, Q\ are bounded and uniformly convex, dilo,d£li £ C 2,1 , p £ 
C 1 ' 1 ^), o- G C 1 ' 1 ^), and A < p(x) < A in Q . Then 

p* G c 3 > a {n ) n c 2 ' a (n ) Voe(o,i), 

anci i/iere exists a constant C which depends only on a, 0,q, fij, Ao, Ai, ||p||<ji,i, 1 1 c" 1 1 c 1 ' 1 suc ^ 
ll-P*llc8,a (no ) < C and ||i"||oa,a(n ) < C. 

Moreover, there exist positive constants c\ and c 2 and k, depending only on Ao, Ai, ||p||co,a, 
and \\a\\(jo,a, such that 

did <X7 2 P*(x) < c 2 Id Viefio 

and 

The first statement is standard optimal transport theory, see |101 [22] , except for the fact that we are 
not assuming that the second moment of \x is finite, thus the classical Wasserstein distance from [i and 
v can be infinite. Nevertheless the existence of an "optimal" map is provided by [22]. The W 2,1 part 
of the second statement follows from a recent regularity result about solutions of the Monge- Ampere 
equation [15], while the C 1 '^ regularity was proven by Cafiarelli in [6j[U[T0]. The regularity up to the 
boundary and the oblique derivative condition of the third statement have been proven by Caffarelli 
[7] and Urbas [25]. 

Remark 2.2. By compactness and a standard contradiction argument, the constants C\ and C 2 in the 
statement (ii) of the previous theorem remain uniformly bounded if S7i varies in a compact class (with 
respect, for instance, to the Hausdorff distance) of convex sets. In particular, let Q,™ be a sequence 
of open convex sets which converges to Oi with respect to the Hausdorff distance and a n a sequence 
of densities supported on fi" with Ai < a n < A\ on which converge to a in L 1 (M 3 ). Then the 
estimates in Theorem I2.1l fii) hold true with constants independent of n. 

Remark 2.3. As already mentioned in the introduction, in statement (ii) the optimal regularity is 
that for every Q (g J7 an d < A < p(x) < A < oo in fi, there exist 7(A, A, Ai,Ai) > 1 and 
C(f2, A, A, Ai, Ai), such that 

I < C. 

Jn 

However, as explained in Remark 13. 6[ this improvement does not give any advantage. 

3. The dual problem and the regularity of the velocity field 

In this section we recall some properties of solutions of (jl.4p . and we show the L 1 integrability of 
the velocity field u t defined in (jl.5p . 

We have the following result whose proof follows adapting the argument of [13] , where compactly 
supported initial data are considered. Since the velocity Ut has at most linear growth, the speed of 
propagation is locally finite and the proof readily extends to general probability densities. 
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Theorem 3.1 (Existence of solutions of (|1.4[> ). Let Pq : K 3 — > R be a convex function such that 
(VPo)jj^h < Then there exist convex functions P t ,P t * : R 3 -> R suc/i t/icrf (VPt) tt J% = Pt^ 3 , 

(VPf)^pt = J^n, Ut{x) = J(x — VP t *(i)), and p t is a distributional solution to (|1.4j) . namely 

(3.1) / / <^ d t (ft{x) + Vift(x) ■ Ut(x) \pt{x) dxdt + / <£o(x)poO*0 = 

/or every v? G C^°(M 3 x [0,oo)). 

Moreover, the following regularity properties hold: 

(i) piJz? 3 G C([0,oo),P w (R 3 )), w/iere V W (M. 3 ) is the space of probability measures endowed with 
the weak topology induced by the duality with Co(R 3 ); 

(ii) P t * - P*(0) G L&([0, oo), w£ c °°(K 3 )) n C([0, oo), ^(R 3 )) /or every r G [1, oo); 

(iii) |C/((x)| < |x| + dn for almost every x G R 3 , /or a// i > 0. 

Observe that, by Theorem 13 .1( h)- t t— > pt^£ 3 is weakly continuous, so pt is a well-defined function 
for every t > 0. Further regularity properties of Pt and P t * with respect to time will be proven in 
Proposition 13.51 

In the proof of Theorem 1 1.31 we will need to test with functions which are merely W ' with compact 
support. This is made possible by a simple approximation argument which we leave to the reader, see 
[TJ Lemma 3.2]. 

Lemma 3.2. Let pt and Pt be as in Theorem \3.1[ Then (|3.1|) holds for every if G M /1 ' 1 (IR 3 x [0, oo)) 
which is compactly supported in time and space, where now (fo(x) has to be understood in the sense of 
traces. 

Lemma 3.3 (Space-time regularity of transport). Let S!C1 3 be a uniformly convex bounded domain 
with dn G C 2 '\ let R>0, and consider p G C°°(P(0,P) x [0,oo)) and U G C™(B(0,R) x [0,oo);R 3 ) 
satisfying 

dtfh + V ■ (U t pt) = mP(0,P) x [0,oo). 
Assume that J B ^ Q R ^podx = 1, and that for every T > there exist At and At such that 

< A T < pt{x) < At < oo V (x, t) G B(0,R) x [0,T]. 

Consider the convex conjugate maps Pt and P 4 * such that (VPt)flJ2fo = pt and (VP/Wt = J^n (unique 
up to additive constants in Q and B(0,R) respectively). Then: 



(i) P t * - f B(0}R) P t * G Li Ploc ([0,oo);C 2 '«(P(0,P))). 

(ii) The following linearized Monge- Ampere equation holds for every t G [0, oo): 



(3.2) 



V • (ptiVPn^dtVP?) = -V • (p t Ut) m B(0, R) 

pti^P^dtVPt -v = on dB(0, R). 



Proof. Observe that because pt solves a continuity equation with a smooth compactly supported vector 
field, j B ^ Q pt dx = 1 for all t. 

Let us fix T > 0. From the regularity theory for the Monge- Ampere equation (Theorem 12.11 applied 
to P t and P*) we obtain that P t G C 3 > a (n) n C 2 ' a (Tl) and P* G C 3 > Q (P(0, R)) n C 2 ' a (B(0,R)) for 
every a G (0, 1), uniformly for t G [0, T], and there exist constants ci, c% > such that 

(3.3) Cl Ld<V 2 P t *(x) < c 2 Ld V(x,t) G P(0,P) x [0,T]. 
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Let h E C 2,1 (M 3 ) be a convex function such that Q = {y : h{y) < 0} and |V/i(y)| = 1 on dfl, so that 
Vh(y) = un(y). Since VP t * E C 1,a (B(0, R)), it is a diffeomorphism onto its image, we have 

(3.4) h{VP t *(x)) =0 V (x, t) E 95(0, R) x [0, T]. 

To prove (i) we need to investigate the time regularity of P^ — j~ B ^ mP*- 

Possibly adding a time dependent constant to Pt, we can assume without loss of generality that 
Ib(o R) = ® f° r an ^ * ne con dition (VP^)^pt = JSfn we get that for any < s,t < T and 
x E 5(0, i?) it holds 

p s (x) - Pt {x) _ det(V 2 P;(x)) - det(V 2 P t *(x)) 
s — t s — t 

i- 1 Sdet, , „ ^a^, « , \ fyJ?(s)-fyi?(s) 



Moreover, from (|3.4p we obtain that on dP>(0, R) 

o = h(vp:(x))-h(yp*(x)) 



s - t 



[ V1i{tVP*(x) + (1 - r)VP t *(x)) dr ■ 
Jo 



s - t 

Now, given a matrix A = (&,•), we denote by M{A) the cofactor matrix of A. We recall that 
(3.T) ^ = A % (A), 

and if A is invertible then M(^4) satisfies the identity 

(3.8) M(A) = det(A)A~ l . 

Moreover, if A is symmetric and satisfies c\ld < A < c%ld for some positive constants c\, C2, then 

(3.9) c\ld < M(A) < c 2 2 Id. 
Hence, from (|3"3]h (13771) . (1331) and ([379]) it follows that 

(3.10) ^—^ = yj f 1 MlJ (rV 2 P: + (1 - r)V 2 P t *) dr) d t ( K " ^ 
with 



c 2 /d < f A%(rV 2 P; + (1 - r)V 2 P;) dr < cpd. 
Jo 



Also, from Theorem I2.1f iii) the oblique derivative condition holds, namely there exists k > such 
that 

Vh(VP t *(x)) • u B{0jR) (x) > K Vi E &B(0,i?). 

Thus, since 



uniformly in f and x, we have that 



lim [ Vh(TVP*(x) + (1 - t)VP?(x)) dr = Vh(VP?(x)) 
s ->tJo 

e have that 

j\h(rVP:(x) + (1 - r)VP t *(x))dr ■ v B(m (x) > \ 
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for s — t small enough. 

Hence, from the regularity theory for the oblique derivative problem |191 Theorem 6.30] we ob- 
tain that for any a S (0,1) there exists a constant C depending only on f2, T, a, \\(p s — pt)/(s — 
t )llc°."(i?(o, J R))> such that 

i?(x)-i?(s) 



s - t 



< C. 



C 2 > a (B(0,R)) 

Since dtpt £ L°°([Q, T], C 0,a (B(0, R))), this proves point (i) in the statement. To prove the second 
part, we let s — > t in (|3.10|) to obtain 

3 

(3.11) d tPt = ^ M i3 (V 2 P;{x)) OtdijP^x). 

Taking into account the continuity equation and the well-known divergence-free property of the co- 
factor matrix 

3 

^ / d i M ij (V 2 P t *(x))=0, j = 1,2,3, 

8=1 

we can rewrite (|3.1ip as 

3 

-v • (u tPt ) = HM^PZix^dtdjPZix)). 

»J=1 

Hence, using (|3.8p and the Monge-Ampere equation det(V 2 i-f) = pt, we get equation (|3.2p . 

In order to obtain the boundary condition in (|3.2p . we take to the limit as s — > t in (|3.6|) to get 

(3.12) X7h(X7P t *(x)) • d t VP t *(x) = 0. 

Since h satisfies 0, = {y : h(y) < 0} and V-Ff maps B(0,R) in Q, we have that B(0,R) = {y : 
h(X7P t *(y)) < 0}. Hence v B (o,R)(x) is proportional to V[h o VP t *](x) = V 2 P t * (x)V h(V P t * (x)) , which 
implies that the exterior normal to O, at point VP t *(s), which is V/i(V?(*(x)), is collinear with 
Pt(V 2 P t *) _1 ^B(o,.R)- Hence from (13TT2"]) it follows that 

^(v 2 ^*)- 1 ^,*) • d t vp; = o, 

as desired. □ 

Lemma 3.4 (Decay estimates on pt). Let vt '■ x [0,oo) — > R 3 be a C°° velocity field and suppose 
that 

sup|V -v t {x)\ < N, \v t {x)\ < A\x\ + D V (x, t) £ R 3 x [0, oo) 

x,t 

for suitable constants N, A, D. Let po ^ e a probability density, and let p t be the solution of the 
continuity equation 

(3.13) d t pt + V • (v t pt) = in R 3 x (0, oo) 
starting from p$ . Then: 

(i) For every r > and t € [0, oo) it holds 

(3-14) \\Pt\\oo < e Nt \\po\\ooi 

At i 

(3.15) p t (x) > e~ m inf |po(y) : y e B (o, re M + D & ~ ) | Vx£5(0,r). 
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(ii) Let us assume that there exist do S [0, oo) and M S [0, oo) such that 

(3.16) Po(x) < j — whenever \x\ > M. 

Then for every t G [0, oo) we have that 

d 2 K e^ N+AK ^ e At - 1 

(3.17) p t {x) < — whenever \x\ > 2Me At + 2D — . 

\x\ n A 



(iii) Let us assume that there exists R > such that po is smooth in B(0,R), vanishes outside 
B(0,R), and that vt is compactly supported inside B(0,R) for all t > 0. Then pt is smooth 
inside B(0,R) and vanishes outside 5(0, i?) for all t > 0. Moreover if 0<X<po<A<oo 
inside B(0, R), then 

(3.18) \e~ tN <p t < Ae tN inside B(0, R) for all t > 0. 

Proof. Let Xt(x) £ C°°(]R 3 x [0, oo)) be the flow associated to the velocity field vt, namely the solution 
to 



(3.19) 



'f t X t (x) = v t (X t (x)) 
X (x) = x. 



For every t > the map 1 1— > Xt(x) is invertible in K 3 , with inverse denoted by Xf 1 . 

The solution to the continuity equation (|3. 13[) is given by pt = X t $po, and from the well-known 
theory of characteristics it can be written explicitly using the flow: 

(3.20) Pt (x) = po(X- 1 (x))eti y - v °( x °( x t 1 W> ds V (x,t) £l 3 x [0,oo). 
Since the divergence is bounded, we therefore obtain 

(3.21) p (X t -\x))e- Nt < p t (x) < p Q {Xt\x))e m 

Now we deduce the statements of the lemma from the properties of the flow Xf. 
(i) From (|3.21|) we have that 

Pt{x) < e Nt p (X- 1 (x)) < e Nt sup Po (x), 

which proves (|3.14p . From the equation (|3.19p we obtain 

£\X t (x)\ < \d t X t (x)\ < A\X t (x)\ + D 
dt 

which can be rewritten as 

(3.22) - A\X t (x)\ — D < j t \X t (x)\ < A\X t (x)\ + D. 
From the first inequality we get 

1 - e~ At 



\X t (x)\ > \x\e~ At - D 



A 

which implies 

At _ i 

\ x \e At + D e -—±>\X T \x) 
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or equivalently 

(3.23) Xf 1 ({\x\ < r}) C < re + D }. 

Hence from (|3.21[) and (|3.23p we obtain that, for every x G P(0,r), 

Pt(x) > e- m p (Xr\x)) 

> e- Nt mf{p o (y):yeXr 1 (B r (0))} 



> e- Nt mf 



{po(y):\y\<re At + D— — } 



which proves (|3.15|) . 

(ii) From the second inequality in (I3.22p . we infer 

e At - 1 

\X t (x)\<\x\e At + D— — , 

which implies 

At _ 1 

(3.24) |z| < \Xt l {x)\e At + D - . 

Thus, if \x\ > 2Me At + 2D^^, we easily deduce from (^24^ that \X^ l (x)\ >M + \x\e~ At /2, so by 

^ C Nt J ryK „(N+AK)t 

Pt{x)<e p (X t (g)) < < ^p? > 

which proves ()3.17p . 

(hi) If = in a neighborhood of <9P(0, R) it can be easily verified that the flow maps Xt : M 3 — >• M 3 
leave both P(0, R) and its complement invariant. Moreover the smoothness of vt implies that also Xt 
is smooth. Therefore all the properties of pt follow directly from (|3.20p . □ 



We are now ready to prove the regularity of VP t * . 

Proposition 3.5 (Time regularity of optimal maps). Let $7 C M 3 be a bounded, convex, open set and 
let du be such that Q C B(0,d$i). Let pt and Pt be as in Theorem \3.1\ in addition let us assume that 
there exist K > 4, M > and cq > such that 

(3.25) po(x) < j — whenever \x\ > M . 

Then VP 4 * G W£ c (M 3 x [0,oo);M 3 ). Moreover for every k G N and T > i/iere exists a constant 

C = C(k,T,M,CQ, \\po\\oo,dn) such that, for almost every t G [0,T] it ZioWs 

(3.26) 

/ p t |^vp;|iogt(|^vp;i)dx<2 3 ^ 1 ) / ft |v 2 p*|iog 2 + fc (|v 2 p;i)^ + c vr>o. 

ifi(0,r) JB(0,r) 
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Proof. Step 1: The smooth case. In the first part of the proof we assume that O is a convex smooth 
domain, and, besides (|3.25j) . that for some R > the following additional properties hold: 



(3.27) Pt G C°°(S(0,i?) x R), U t G C™(B(0,R) x M;IR 3 ), |V ■ U t \ < N 

(3.28) Al B(0)jR) (x) < < Al B (o,fl)(ic) Vx G M 3 , 

(3.29) + V • (U tP t) = inR 3 x [0, oo) , 

(3.30) (VP;) ( ft=^, 

(3.31) |t^(x)| < |x| +d n 

for some constants iV, A, A, and we prove that (|3.26|) holds for every t G [0, T\. Notice that in this step 
we do not assume any coupling between the velocity Ut and the transport map VP/ 1 . In the second 
step we prove the general case through an approximation argument. 

Let us assume that the regularity assumptions (|3.27p through (|3,3ip hold. By Lemma 13.41 we infer 
that, for any T > 0, there exist positive constants Ay, Ay, ct, Mt, with Mt > 1, such that 

(3.32) *t1b(o,r)(v) < Pt(x) < A t 1b(o,r)(x), 

(3.33) p t {x) < for \x\ > M T , for all t G [0,T]. 



By Lemma [3731 we have that d t Pt G C 2 (B(0, R)), and it solves 

'V • (ptiVP^dtVPt*) = -V • (p t U t ) in B(0,R) 

Multiplying (|3.34p by dtP* and integrating by parts, we get 



(3.34) 



pti^P^dtVPt* -u = in dB(0, R). 



(3.35) 



I p t \(v 2 p t *r 1 / 2 d t vp t *\ 2 dx= [ ptdtvp; -{v 2 p:y l d t vp t 

Jb(q,r) Jb(o,r) 

= - [ ptdtVP t * -U t dx. 
Jb(o.r) 



t dx 



(Notice that, thanks to the boundary condition in (|3.34p . we do not have any boundary term in (|3.35p .) 
From Cauchy-Schwartz inequality it follows that the right-hand side of f)3.35j) can be rewritten and 
estimated by 



f PtdtVP; ■ (V 2 P t *)~ 1/2 (V 2 P t *) 1/2 U t dx 

JB(0,R) 

(3.36) / v 1/2 / x 1/2 

< If Pt \{v 2 p;r l i 2 d t vp;\ 2 dx) If Pt \(v 2 pn l/2 u t \ 2 dx 

\Jb(o,r) J \Jb(o,R) j 

Moreover, the second term in the right-hand side of (|3.36p is controlled by 

(3.37) f p t U t -V 2 P t *U t dx < max (pl /2 \Ut\ 2 ) [ p] /2 \V 2 P t *\ dx. 

JB{Q,R) B(fl,R) V / Jb^R) 

Hence from (|3.35p . (|3.36p . and (|3.37|) we obtain 

(3.38) f p t \(V 2 P t *y 1/2 d t VP t *\ 2 dx < max (p] /2 \U t \ 2 ) f p] /2 \V 2 P t *\dx. 
Jb(o,r) b(o,r) V / Jb(o,r) 



12 L. AMBROSIO, M. COLOMBO, G. DE PHILIPPIS, AND A. FIGALLI 

^From (13.31 p . (I3.32j) . and (|3.33p we estimate the first factor as follows: 

(3.39) max ( p]' 2 (x)\U t {x)\ 2 ) < k)l 2 {M T + d n ) 2 , 

\x\<M T V / 

(3.40) max U ,2 {x)\U t {x)f) < max J t^(\x\ + dn) 2 \ , 

M T <\x\ V / M T <\x\ I 

and the latter term is finite because Mt > 1 and K > 4. 

In order to estimate the second factor, we observe that since V 2 P 4 * is a nonnegative matrix the 
estimate \V 2 P t *\ < AP t * holds (here we are using the operator norm on matrices). Hence, by ()3.32|) 
and ()3.33p we obtain 



AP t *dx ds, 



p] /2 \V 2 P t *\dx < / p} /2 \V 2 P*\dx + p] /2 \V 2 P t *\dx 

B(0,R) J{\x\<M T ] J{\x\>M T } 

< [ A^ 2 AP*dx+ [ -^prAP*dx. 

J{\x\<M T } J{\x\>AI T } \x\ K ' 2 

The second integral can be rewritten as 

/*oo r 

Jo J{\x\>M T }n{\x\- K / 2 >s} 

which is bounded by 

»[A/ T ]- A ' /2 r 

ds / AP^dx. 

J{\x\<s~ 2 / K } 

From the divergence formula, since |VP f *(a;)| < dn (because \7P^(x) G $7 for every x G K 3 ) and 
M T > 1 (so \M T Y K I 2 < 1) we obtain 

(3.41) 

r r r[M T ]- K / 2 r 

/ pl /2 \V 2 P t *\ dx < h^ 2 / \VP t *\d7i 2 + Vcr / ds | VP/ 1 d^ 2 

Jb(o,r) J{\x\=m t } Jo J{\x\=s- 2 / k } 

<ATrA}/ 2 M^dn + ^^dn [ s~ 4/K ds 

Jo 

for all t G [0, T]. Since K > 4 the last integral is finite, so the right-hand side is bounded and we 
obtain a global-in-space estimate on the left-hand side. 

Thus, from (pT38l) . (pT39l) . (|51Sj) . and ([5^I]> . it follows that there exists a constant C x = C\{T, M, c , A, 
(notice that the constant does not depend on the lower bound on the density) such that 

(3.42) / p t \(V 2 P t *y^ 2 d t VP t *\ 2 dx < Ci. 

Applying now the inequality 



(3.43) ablog k + (ab) < 2 k ~ l 



+ 1 

e 



b 2 + 2 3 ( fc - V log? (a) V (a, 6) G R + x 
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(see m Lemma 3.4]) with a = ((V 2 ^*) 1 / 2 ] and b = \(V 2 P?)~ l l 2 d t VP t *{x)\ we deduce the existence of 
a constant C2 = C2(k) such that 

|d t vp t *|iog^(|d t vp t *|) < 2 3 ( fc - 1 )|(v 2 p;) 1/2 | 2 io g 2 f fc (|(v 2 p;) 1/2 | 2 ) + c 2 |(v 2 p;)~ 1/2 9 t vp;i 2 

= 2 3 ( fc ~ 1 )|v 2 p;i io g 2fc (|v 2 p t *|) + c 2 \(v 2 p t *y 1/2 d t vp t *\ 2 . 

Integrating the above inequality over B(0,r) and using (]3.42p . we finally obtain 



B(0,r) 



p t \d t VP t *\log k + (\d t VP*\)dx 



(3.44) <2 3{k ~V f Pt \V 2 P t *\log 2 + k (\V 2 P t *\)dx + C 2 f pt^PZY^dtVPtfdx 

JB(0,r) JB(0,R) 

< 2 3 ( fc ~ 1} f Pt\V 2 P t *\ log^(|V 2 P;|) dx + d • c 2 , 

JB(0,r) 



for all < r < R. 



Step 2: The approximation argument. We now consider the velocity field U given by Theorem 13.11 
we take a sequence of smooth convex domains £l n which converges to £1 in the Hausdorff distance, and 
a sequence (ip n ) C C£°(B(0, n)) of cut off functions such that < i/) n < 1, ip n (x) = 1 inside B(0,n/2), 
\Vip n \ < 2/n in ]R 3 . Let us also consider a sequence of space-time mollifiers o~ n with support contained 
in B(0, 1/n) and a sequence of space mollifiers (p n . We extend the function Ut for t < by setting 
C/j = for every t < 0. 

Let us consider a compactly supported space regularization of po and a space-time regularization 
of J7, namely 

(p * V n ) 



Po 



-l B (0,n), C/f(x) :=([/* (7^, 



where c n t 1 is chosen so that p$ is a probability measure on M 3 . Let p™ be the solution of the continuity 
equation 

d t p n t + V • (U^p?) = in M 3 x [0, 00) 

with initial datum pg . Prom the regularity of the velocity field and of the initial datum p^ we have 
that p n £ C°°(B(0,n) x [0,oo)). 

Since Ut is divergence- free and satisfies the inequality |i7((x)| < |a;| + dn, we get 

\U?\(x) <\U* a n \(x) < \\U t \\ L °o( BMn)) < \x\ + d n + - < \x\ +d Q + l, 

|v • ur\(x) = \(u t * a 11 ) ■ w> n |(x) < 2 ( n + 1 +^) < 3 

n 

for n large enough. Moreover, from the properties of po we obtain that, for n large enough, 

2c ^ 4c 



\Po\ 



< 2||po |oo 



and 



A" 


< | 




1 




Po 7 



X 



K 



V|s| > M + 2, 



< 



L°°(B(0,n)) 



1 

PO 



L°°(B(0,n+l)) 



Hence the hypotheses of Lemma l3.4l are satisfied with N = 3, A = 1, D = c2q+1, do = 4co- Moreover 

-1 



p™ vanishes outside B(0,n), and by (|3.18p there exist constants A n := e 



-3T 



L°°(B(0,ra+l)) 



> 0, 
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A := 2e 3T ||/9o||oo> and Mi, c\ depending on T, M,co,d,Q only, such that 

A„,t < Pt{x) < A V (x, t) G 5(0, n) x [0, T], 

Pt( x ) < j — whenever |x| > Mi. 

(Observe that A„ depends on n, but the other constants are all independent of n.) Thus, from 
Statement (ii) of Lemma 13.41 we get that, for all r > 0, 

(3.45) pUx) >e- 3T inf{^(y) : y G b(o, re* + (do + l)[e* - 1]) } V(M) G 5(0, r) x [0,T]. 

If n is large enough, the right-hand side of f|3.45j) is different from 0, and can be estimated from below 
in terms of po by 

A = A(r, T, po, fi) := e" 3T inf {p (y) : y G B (o, re* + (d n + 1) [e* - 1] + l) } > 0. 

Therefore, for any r > we can bound the density p n from below inside 5(0, r) with a constant 
independent of n: 

(3.46) A < ft(x) < A V (x, t) G 5(0, r) x [0, T]. 

Let now 5 t n * be the unique convex function such that P t n *(0) = and (VP™)jj/9™ = From the 

stability of solutions to the continuity equation with BV velocity field, (3j Theorem 6.6], we infer that 

(3.47) ft -> p t in ^(K 3 ), for any t > 0, 

where pt is the unique solution of (|1.4p corresponding to the velocity field U. S ince r^ n is converging 
to f2, from standard stability results for optimal transport maps (see for instance \26\ Corollary 5.23] 
and [16l Section 4]) it follows that 

(3.48) VPr -> VP t * in Ll c (R 3 ) 

for any t > 0. Moreover, by Theorem 12. lf ii). Remark I2.2[ and (j3.46|) . for every k G N 

(3.49) limsup / p^|V 2 P i ; n *|log^(|V 2 P i ; n *|)(ix< oo Vr > 0, 

rwoo JB(0,r) 

and by the stability theorem in the Sobolev topology estabilished in [161 Theorem 1.3] it follows that 

(3.50) lim / ^|V 2 Pr|log 2 , fc (|V 2 Pri)^= / p t \V 2 P t *\logl k (\V 2 P t *\)dx Vr > 0. 
n ->°° JB(0,r) JB(p,r) 

Since (pt,U^) satisfy the assumptions (|3.27p through (|3.31j) . by Step 1 we can apply f|3.44j) to 
(pj?, U?) to obtain 

(3.51) f p?\d t VP t n *\log k + {\d t VP t n *\)dx < 2 3 ( fc " 1 ) f Pt \V 2 P t n *\logf(\V 2 Pr\)dx + C 

JB(0,r) JB(0,r) 

for all r < n, where the constant C does not depend on n. 

Let (j) G C£°((0, T)) be a nonnegative function. From the Dunford-Pettis Theorem, taking into 
account (pHT)) and (|515J) . it is clear that 4>{t) p? d t V P?* converge weakly in L 1 (P(0,r) x (0,T)) to 
(f>(t) ptdtV P£ . Moreover, since the function k; h- > \w\ log + (\w\/r) is convex for every r G (0, oo), we can 
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apply Ioffe lower semicontinuity theorem [21 Theorem 5.8] to the functions cf){t) p™ dtV P™* and 4>(t)p r t 

to infer 

(3.52) 

f 4>{t) f Pt \d t VP?\iog k + (\d t vp t *\)dxdt<]imM f <p(t) [ P ?\d t vpr\iog k + (\d t vpr\)dxdt. 

JB(0,r) n ^°° J0 JB(0,r) 

Taking (|3.5ip . ()3.50p . (|3.49p . and (|3.52p into account, by Lebesgue dominated convergence theorem 



we obtain 



(f>(t) [ Pt \d t vp t *\iog k + (\d t vp t *\)dxdt 

JB(0,r) 

[ T <f>(t) [2 3 ^ [ P t\v 2 p;\ io g ^(|v 2 p;i) dx + c ) dt. 

JO \ JB(0,r) 



JB(0,r) 

< 



Since this holds for every G C^°((0, T)) nonnegative, by a localization argument we obtain the 
desired result. □ 

Remark 3.6. Thanks to Remark 12.31 one can prove that for every T > and r > there exist a 
constant k > 1 and a constant C which depend on r, po, T, dn such that, for almost every t £ [0,T] 
we have that VP t * £ W 1 ' K (B(0,r) x [0,oo);IR 3 ) and 

/ p t \d t VP t *\ K dx < C. 

JB(0,r) 

This estimate provides better local integrability of the time derivative of VP*. The proof follows the 
same lines of Proposition 13.51 (see also [201 Proposition 5.1]). However the exponent k is not universal, 
but depends in a nontrivial way from the local lower bounds on the density which are related to r, po 
and T. Therefore we preferred to state Proposition 13.51 with a universal modulus of integrability. 

We finally point out that in the compact setting studied in [1] the same argument provides a global 
L K estimate of dt^Pt on the torus, with k depending only on the upper and lower bound on po, which 
is also uniform in time. 

4. Existence of an Eulerian solution 

Proof of Theorem First of all notice that by statement (ii) of Theorem 12.11 and Proposition 13.51 
it holds \V 2 P t *\, \d t VP t *\ G L^ c ([0, oo), Ll oc (R 3 )). Moreover, since (VP t ) tt Jz? 3 = p t , it is immediate to 
check the function u in (fT3D is well-defined and \u\ belongs to L£ c ([0, oo), Lj^QR 3 )). 

Let 4> £ C£°(Q x [0, oo)) be a test function and let us consider ip : M 3 x [0, oo) — > M 3 given by 

(4.i) Mv) ■= vMVPtiv))- 

Clearly ip is compactly supported in time because so is <f>; moreover P t are Lipschitz on supp 4>t as t 
varies in any compact subset of [0, oo) with bounded Lipschitz constants. Hence the set V Pt(supp (j>t) , 
which contains suppy^, is bounded in space. Therefore 4>t is compactly supported in IR 3 x [0, oo). 
Moreover, Proposition 13.51 implies that (p E VF 1 ' 1 (M 3 x [0, oo)). So, by Lemma 13.21 each component of 
the function <pt(y) is an admissible test function for (|3.ip . For later use, we write down explicitly the 
derivatives of <p: 



(4.2) 



'd m (y) = y[d t <j> t ](VP t *(y)) + y([V&](i?(y)) • d t VP t *(y)), 
V<pt(y) = ldMVP t *(y)) + y®{lV T &](P t *(y))V 2 P t *(y)). 
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Taking into account that (VPi)ji^ = pt^f 3 and that [VP/](VPt(x)) = x almost everywhere, we can 
rewrite the boundary term in (|3.ip as 



(4.3) / My)Mv)dy = / VP (x)Mx)dx. 

Jr 3 Jn 

In the same way, since Ut(y) = J(y — VP t *(i/)), we can use (|4.2|) to rewrite the other term as 
{d m (y)+V<p t (y)-U t (y)} 

= I™ I {^Pt(x)dtU^) + ^Pt(x){V(t>t(x)-[d t VP;)(VP t (x))) 
Jo Jn 1 

+ [id<k(x) + vp t (x) ® (v T Mx)v 2 p t *(yp t (x)))]j(yp t (x) - s)} 



(4.4) 



<ix (it 



which, taking into account the formula (jl.5p for u, after rearranging the terms turns out to be equal 
to 

(4.5) J JvP t (x){dtMx)+Mx)-VMx)} + j{vPt(x)-x}<t> t (x)dxdt. 

Hence, combining ([P]) . (|4T5|) . and ([3TT]1 . we obtain the validity of (fL6j) . 

Now we prove ([T7f|l . Given <E C c °°(0,oo) and ip £ C^°(0), let us consider 93 : R 3 x [0, 00) ->■ R 
defined by 

(4.6) ^(y) :=0(i)^(VP t *(y)). 

As in the previous case, </? € W 1,1 (R 3 x [0, 00)) and is compactly supported in time and space, so we 
can use ip as a test function in (|3.ip . Then, identities analogous to (|4.2j) yield 



+ 



= / / {d m (y) + V<p t (y)-U t (y)}pt(y)dydt 

/ ip{x)dxdt 
Jn 

J <f>(t) j {v^(x) • d t VP t *(VP t (x)) + V T ^j(x)V 2 P t *(VP t (x))J(VP t (x) - x)} tfccft 

0(t) / Vip(x) • u t {x)dxdt. 
Jn 

Since (ft 1S arbitrary we obtain 

/ Vip(x) ■ ut(x) dx = for a.e. t > 0. 
in 

By a standard density argument it follows that the above equation holds outside a negligible set of 
times independent of the test function ip, thus proving (|1.7|) . □ 
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